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Abstract 

Histone methylation plays an important role In gene transcription and chromatin organization and Is linked to the silencing 
of a number of critical tumor suppressor genes in tumorlgenesis. G9a is a histone methyltransferase (HMTase) for histone H3 
lysine 9. In this study, we investigated the role of G9a In neuroblastoma tumor growth together with the G9a Inhibitor 
BIX01294. The exposure of neuroblastoma cells to BIX01294 resulted In the Inhibition of cell growth and proliferation, and 
BIX01294 treatment resulted In the Inhibition of the tumorigeniclty of neuroblastoma cells In NOD/SCID mice. Therefore, 
G9a may be a potential therapeutic target in neuroblastoma. Moreover, we found several specific characteristics of 
autophagy after BIX01294 treatment, including the appearance of membranous vacuoles and microtubule-associated 
protein light chain 3 (LC3B). Similar results were observed In G9a-knockdown cells. In conclusion, our results demonstrated 
that G9a is a prognostic marker in neuroblastoma, and revealed a potential role of G9a In regulating the autophagy 
signaling pathway In neuroblastoma. 
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Introduction 

Tumorigenesis is considered to be a multi-step process ranging 
from stages characterized by normal histological features to 
carcinoma features. Epigenetics has been recently defined as 
inheritable changes in gene expression not due to any alteration in 
the DNA sequence. Histone methylation is the fundamental 
epigenetic mechanism that regulates gene expression in cancers 
and is linked to the silencing of a number of critical tumor 
suppressor genes in tumorigenesis [1,2]. Recently, G9a was 
reported to be a major H3K9mel and H3K9me2 HMT in vivo 
[3-6] , and several studies have identified the critical role that G9a 
plays in various biological processes, including embryo develop- 
ment, immune response, drug response and tumor cell growth [7- 
14]. Moreover, current evidence suggests that G9a promotes 
invasion and metastasis in lung cancer [13], and highly expressed 
G9a was observed in hepatocellular carcinomas [15]. Therefore, 
G9a may be a key regulator that serves as a potential therapeutic 
target during tumor formation. 

In addition, autophagy is an evolutionarUy conserved mecha- 
nism that involves the degradation of macromolecules, ribosomes, 
and organelles [16]. Autophagy is the primary intracellular 
catabolic process responsible for long-lived protein and organelle 
degradation and recycling, whereas the ubiquitin/proteasome 
system is the major cellular pathway responsible for short-lived 



protein degradation [17,18]. The following four primary forms of 
autophagy have been described: macroautophagy (referred to here 
as autophagy), selective autophagy, microautophagy, and chaper- 
one-mediated autophagy [19-21]. Autophagy serves as an 
adaptive response to cellular stress such as hypoxia and nutrient 
deprivation, which involves the synthesis of a double-membrane 
structure known as the phagophore. The phagophore ultimately 
elongates and closes to sequester cytoplasmic proteins and 
organelles, forming the autophagosome, and undergoes a stepwise 
maturation process [22-24]. Mammalian autophagy-related genes 
(ATG) participate in distinct steps of autophagy. For example, 
microtubule-associated protein light chain 3 (LC3B) undergoes 
lipidation and is recruited to the phagophore where it is essential 
for membrane elongation and closure [20]. 

Neuroblastoma is a common childhood malignant tumor of 
neural crest origin, arising in the sympathetic nervous system, and 
this condition accounts for approximately 10% of pediatric 
cancers and 15% of cancer-related deaths in children [25-31]. 
In this study, we investigated the role of G9a in neuroblastoma 
tumor growth together with BIX01294, which is a specific G9a 
inhibitor [32-34]. We provide experimental evidence supporting 
the role of G9a in the transcriptional regulation of autophagy in 
neuroblastoma cells. Pharmacological inhibition or RNA interfer- 
ence (RNAi) of G9a led to increased LC3B expression and 
autophagosome formation. Collectively, we identified G9a as a 
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Figure 1. High G9a expression is a prognostic indicator of poor survival in neuroblastoma patients. A, Kaplan-Meier analysis of 
progression-free survival for the Versteeg database with the log rank test P value indicated. B, box plot of G9a expression levels in stage (ST) 1-4S 
tumors. C, box plot of G9a expression levels in tumors from patients in the non-death and tumor-caused-death groups. D, box plot of G9a expression 
levels in tumors from patients in the <18-month and >18-month groups. E, Kaplan-iVleier analysis of progression-free survival for the Neuroblastoma 
Prognosis Database w\th the log rank test P value indicated. F, box plot of G9a expression levels in tumors from dead and alive groups. G, western 
blot analysis of G9a expression in five neuroblastoma cell lines. For the data in A, the G9a cutoff value of 7.62 was used to separate the patients into 
high and low G9a expression groups. The data in B (ST1 vs. ST2, ST2 vs. ST3, ST3 vs. ST4, and ST4 vs. ST4S), C, D, and F were analyzed using two-tailed 
student's t-test with the P values indicated. For the data in E, the G9a cutoff value of 0 was used to separate the patients into high and low G9a 
expression groups. 

doi:1 0.1 371/journal.pone.01 06962.g001 
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Figure 2. Inhibition of G9a represses neuroblastoma cell growth and proliferation. A, morphologic examination of five neuroblastoma 
cells treated with BIX01294 or water for 2 and 4 days. Scale bars, 5 |im. B, neuroblastoma cells were either treated with 5 jilVI BIX01294 or water for 2 
days and 4 days, respectively, then analyzed for cell counting with the TC10 Automated Cell Counter, error bars, SD, n = 5. Statistical analysis was 
performed using two-tailed student's t-test, *p£0.01. C, neuroblastoma cells were either treated with 5 |ilVl BIX01294 or water and then analyzed for 
cell growth by the CCK8 assay. Each value represents the average obtained from five independent experiments; error bars, SD. Statistical analysis was 
performed using two-tailed student's t-test, *p£0.01. 
doi:1 0.1 371 /journal.pone.01 06962.g002 



prognostic marker for survival in patients with neuroblastoma and 
a regulator of neuroblastoma cell growth, proliferation and 
autophagy. Our results suggest a novel potential role of G9a in 
the regulation of the autophagy signaling pathway in neuroblas- 
toma. 

Materials and Methods 

Cell culture 

The neuroblastoma cell line BE(2)-C was grown in a 1:1 
mixture of Dulbecco's modified Eagle's medium and Ham's 
nutrient mixture F12 (DMEM/F-12) (Life Technologies, Grand 
Island, NY, United States) supplemented with 10% fetal bovine 
serum (FBS) (Life Technologies). Other neuroblastoma cells (SK- 
N-AS, SK-N-DZ, SK-N-Fl, and SHEPl) were cultured in 
Dulbecco's modified Eagle's medium (DMEM) (Life Technologies) 
plus 10% FBS. AU of the cells were obtained from the American 
Type Culture Collection (ATCC, Manassas, VA, United States) 
and cultured at 37°C in a 5% CO2 humidified incubator. 

Cell proliferation and cell cycle assays 

The G9a inhibitor BIX01294 (B9311, Sigma-Aldrich, St. Louis, 
MO, United States) was dissolved in water. The cells were grown 
to 60-70% confluence and treated with the indicated concentra- 
tions of BIX01294 [35]; cells treated with water were used as a 
control. After treatment, the adherent and floating cells were 
pooled at different time points, collected by centrifugation, and 



then washed once with ice-cold PBS. The sample obtained was 
analyzed with the TCIO Automated Cell Counter (Bio-Rad, 
Hercules, CA, United States), and cell counting was determined 
using trypan blue dye (#145-0021, Bio-Rad). The samples for the 
cell cycle assay were fixed with 70% ethanol, stained with PI, and 
analyzed by flow cytometry (BD FACSVerse, BD Biosciences, 
San Jose, CA, United States). The data were analyzed with 
CellQuest Pro software (BD Biosciences). The cell growth curve 
was detected by CCK-8 (CK04-05, DOJINDO, Kamimashiki 
gun, Kumamoto, Japan). 

Immunofluorescence analysis 

For immunofluorescent staining, the cells were grown on 
coverslips. After treatment with BIX01294, the cells were washed 
with PBS and fixed in 4% paraformaldehyde in PBS for 20 min at 
room temperature. Then, the cells were permeabUized with 0.3% 
Triton X-100 for 5 min. The cells were blocked with 5% milk, 
incubated with a primary antibody, and then incubated with the 
appropriate secondary antibody. The primary antibodies were 
used at a dilution of 1:500 for the rabbit monoclonal antibody 
against LC3B (Cell Signaling Technologies, Danvers, MA, United 
States) and 1:2000 for the mouse monoclonal antibody against ot- 
tubulin (clone B-5-1-2, Sigma-Aldrich). Alexa Fluor 488 Goat 
Anti-Rabbit IgG (H-hL) and Alexa Fluor 594 Goat Anti-Mouse 
IgG (H-hL) (Life Technologies, Carlsbad, CA, United States) were 
used as secondary antibodies. In total, 300 nM DAPI in PBS was 
used for nuclear staining. 
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Figure 3. Inhibition of G9a induces cell cycle arrest in G1 phase. A and B, neuroblastoma cells were either treated with 5 |ilVl B1X01294 or 
water for 2 days and analyzed for the cell cycle by flow cytometry. Each column represents the average obtained from three independent 
experiments; error bars, SD. Statistical analysis was performed using two-tailed student's t-test, *p<0.01 . C, western blot analysis of cyclins and CDKs 
associated with G1 phase in neuroblastoma cells treated with BIX01294 or water for 2 days. a-Tubulin levels are shown as the loading control. 
dol:1 0.1 371/journal.pone.01 06962.g003 



Lentiviral infection 

The lentiviral human G9a shRNA constructions were obtained 
from Open Biosystems (Thermo Fisher Scientific, Pittsburgii, PA, 
United States). The lentiviruses were produced by co-transfection 
with an shRNA-expressing vector or GFP shRNA. The lentiviral 
constructs were transfected into 293FT packaging cells using 
Lipofectamine 2000 (Life Technologies). Virus-containing super- 
natants were harvested and titered and then were used to infect 
target cells with 4 |ig/ml polybrene (sc- 134220, Santa Cruz, 
Dallas, Texas, United States). One day after the final round of 
infection, the cells were cultured in the presence of 2 ng/ml 
puromycin (All 13803, Life Technologies) for 3 days, and the 
drug-resistant cells were pooled. 

Western blot analysis 

Western blot analysis was performed with the primary 
antibodies including anti-G9a (1:500, Santa Cruz), anti-a-tubulin 
(1:2000, Sigma-Aldrich), anti-H3K9me2 (1:,500, Abeam, Cam- 
bridge, MA, United States), anti-CDK2 (1:500, Santa Cruz), anti- 
CDK4 (1:500, Santa Cruz), anti-CDK6 (1:500, Santa Cruz), anti- 
CyclinDl (1:500, Santa Cruz), anti-CyclinE (1:500, Abeam), anti- 
LC3B (1:500, Cell Signaling Technologies), and the autophagy 
antibody sampler kit (1:500, Cell Signaling Technologies). 
Horseradish peroxidase-conjugated goat anti-mouse (1:20000), 
goat anti-rabbit (1:20000), and rabbit anti-goat (1:10000) immu- 
noglobulin G (IgG; KPL, Gaithersburg, Maryland, United States) 
were used as secondary antibodies. 



Soft agar clonogenic assay 

In total, 1 xlO'^-2.5xlO'^ cells were mixed with 0.3% Noble 
agar in growth medium and plated into six-well plates containing a 
solidified bottom layer {0.6% Noble agar in growth medium). The 
colonies were photographed after 14 to 21 days and recorded. 

In vivo tumorigenic assay 

Six age-matched non-obese diabetic severe combined immuno- 
deficient (NOD/SCID) female mice (4 weeks old) were used in this 
study. The mice were housed with 2-3 animals per cage and 
maintained under Specific Pathogen-Free (SPF) conditions. For 
the experimental tumorigenesis assays, 1 x 10^-3x10 cells were 
resuspended in 100 |j.l DMEM and injected subcutaneously into 
both flanks of each mouse. After 1 week of tumor growth, the mice 
were randomly divided into two groups. One group was injected 
intraperitoneally with BIX01294 at 4 mg/kg (mice body weight), 
and the other group was injected with water as a control. Tumor 
growth was measured by caliper measurement, and tumor volume 
was calculated with the formula 4/3nr^, where r is the radius of the 
tumor. Four weeks after injection, the tumors were removed and 
weighed. 

Etiiics statement 

This study was conducted in accordance with the approved 
guidelines. The protocol was pre-approved by the Institutional 
Animal Care and Use Committee of Southwest University. AH 
efforts were made to minimize the suffering of the animals. 
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Figure 4. Inhibition of G9a induces autophagy in neuroblastoma cells. A, immunofluorescence analysis of five neuroblastoma cells treated 
with BIX01 294 or water for 2 and 4 days. Scale bars, 5 |j.m. B, Statistical analysis of LC3B puncta in neuroblastoma cells treated with BIX01 294 or water 
for 2 and 4 days. Each column represents the average obtained from three independent experiments; error bars, SD. Statistical analysis was 
performed using two-tailed student's t-test, *p<0.01 . C, western blot analysis of G9a function in neuroblastoma cells treated with BIX01 294 or water 
for 2 days. D, western blot analysis of autophagy-related genes in neuroblastoma cells treated with BIX01294 or water for 2 days. a-Tubulin levels are 
shown as the loading control. 
doi:1 0.1 371 /journal.pone.01 06962.g004 



Patient data analysis 

Patient data and gene expression datasets were obtained from 
the Oncogenomics Section Data Center (http://pob.abcc.ncifcrf. 
gov/cgi-bin/JK) and R2: microarray analysis and visualization 
platform (http:/ /hgserverl. amc.nl/cgi-bin/r2/main.cgi). AH prog- 
nosis analyses were conducted online, and all data and P values 
(log-rank test) were downloaded. Kaplan-Meier analysis and the 
resulting survival curves were performed using GraphPad Prism 
(version 6.0). All cutoff values for separating high and low 
expression groups were determined by the online R2 or 
Oncogenomics database algorithm [36,37]. 

Statistical analysis 

AU observations were confirmed by at least three independent 
experiments. Quantitative data are expressed as the mean ±SD. 
Two-tailed Student's t-test was performed for paired samples. P< 
0.05 was considered statistically significant. 

Results 

G9a expression in neuroblastoma is associated with poor 
prognosis 

To investigate the possibility of G9a as a prognostic marker in 
neuroblastoma, we conducted a microarray-based search using the 



Tumor Neuroblastoma public - Versteeg database [36], which is 
available from the online R2: microarray analysis and visualization 
platform. The Versteeg database contains a cohort of 88 patients 
with neuroblastoma representative of various tumor stages and 
genetic alterations. Kaplan-Meier analysis of progression-free 
survival for the Versteeg database showed that high G9a 
expression was strongly associated with a poor outcome, whereas 
low G9a expression was correlated with good overall survival 
(Fig. lA). Moreover, G9a expression was different in various 
tumor stages (Fig. IB), and G9a expression increased significantiy 
in stage 4 tumors compared to stages 3 and 4S (Fig. IB). Then, we 
examined the correlation of G9a expression levels with patient 
cause of death using data available from the Versteeg database. 
The result revealed that the G9a expression level was significandy 
higher in tumor-caused death than non-death group (Fig. IC). 

A younger patient at diagnosis typically has a good prognosis 
[38]; therefore, we examined the correlation of G9a expression 
levels with patient age at diagnosis using information available 
from the Versteeg database. We chose the cutoff at 18 months of 
age because it provided the most significant P value for prognostic 
evaluation [39]. The younger age group expressed significandy 
lower levels of G9a than the older age group (Fig. ID), showing a 
correlation of low G9a expression with younger age at diagnosis. 
Overall, our data suggest that higher levels of G9a predict poor 
prognosis in neuroblastoma. We confirmed that the high 
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Figure 5. Inhibition of G9a decreases tKie tumorigenicity of neuroblastoma cells. A, neuroblastoma cells were plated at 1 xlO to 2.5x10 
cells per well in six-well culture plates. After 14 to 21 days of culture, soft agar colonies developed with cells treated with water. As shown, the cells 
treated with 5 |.iM BIX01294 were observed to give rise to small and scanty colonies in soft agar. Scale bars, 50 nm. B, colonies that were larger than 
0.5 mm or that contained more than 50 cells were recorded. Each column represents the average obtained from three independent experiments; 
error bars, SD. Statistical analysis was performed using two-tailed student's t-test, *p<0.01. C, tumor growth in NOD/SCID mice injected with the 
indicated neuroblastoma cells as measured by caliper measurements. D, scatter plot of xenograft tumor weight with horizontal lines indicating the 
mean per group. The data were analyzed with 2-tailed Student t test, and the P value is indicated. 
doi:1 0.1 371/journal.pone.0106962.g005 



expression of G9a is prognostic of unfavorable outc:omes with the 
Neuroblastoma Prognosis Database available from the online 
Oncogenomics database [37], which includes a cohort of 56 
neuroblastoma patients (Fig. IE and F). Together, our analyses of 
two independent microarray databases indicate that G9a is a 
prognostic marker in neuroblastoma. 

Then, we confirmed the G9a expression levels in five 
neuroblastoma cell lines, SK-N-AS, BE(2)-C, SK-N-DZ, SK-N- 
Fl, and SHEPl (Fig. IG). We found that G9a is highly expressed 
in all five cell lines, indicating that G9a is commonly expressed in 
neuroblastoma. 

G9a inhibition represses neuroblastoma cell growth and 
proliferation 

We next examined the functional consequence of high 
expression of G9a in neuroblastoma cells. The SK-N-AS, BE(2)- 
C, SK-N-DZ, SK-N-Fl, and SHEPl neuroblastoma cell lines 
were treated with BIX01294, a specific inhibitor of G9a. As shown 



in Figure SI A, G9a inhibition suppressed BE(2)-C cell prolifera- 
tion in a dose-dependent manner. AU five cell lines were very 
sensitive to BIX01294 treatment, which dramatically repressed cell 
proliferation after 4 days of treatment (Fig. SIB). This result was 
confirmed by the cell counting kit-8 (CCK8) assay (Fig. 2C), and 
the cell number was highly decreased after treatment with 5 |iM 
B1X01294 (Fig. 2B). We further examined the cell cycle in 
neuroblastoma cells. After drug treatment, the number of cells in 
the Gl phase increased and the number of cells in the S phase 
decreased in SK-N-AS, BE(2)-C, SK-N-DZ, SK-N-Fl, and 
SHEPl cultures (Fig. 3A and B), suggesting that G9a inhibition 
caused by BIX01294 treatment induced cell cycle arrest in Gl 
phase. Western blot analysis showed that G9a inhibition led to a 
marked down-regulation of CycfinDl, CDK4 and CDK6, which 
are collectively required for the cell cycle progression from Gl to 
the S phase [40,41]. Moreover, G9a inhibition led to an obvious 
down-regulation of CDK2 and no significant change in CyclinE, 
which are required for the cell cycle progressions from S to G2 
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Figure 6. Downregulation of G9a represses neuroblastoma cell growth and proliferation. A, western blot analysis of G9a in 
neuroblastoma cells BE(2)-C expressing GFPsi or individual G9asi sequences. a-Tubulin levels are shown as the loading control. B, morphological 
examination of three neuroblastoma cells expressing GFPsi, G9asi#3 or G9asi#4, respectively. The cells expressing GFPsi are shown as the biological 
controls. Scale bar, 5 |im. C, neuroblastoma cells expressing GFPsi, G9asi#3 or G9asi#4 were analyzed for cell counting using the TC10 Automated 
Cell Counter, error bars, SD, n = 5. Statistical analysis was performed using two-tailed student's t-test, *p£0.01. D, proliferation assays of three 
neuroblastoma cells expressing GFPsi, G9asi#3 or G9asi#4, error bars, SD, n = 5. E, western blot analysis of cyclins and CDKs related to the G1 phase 
in neuroblastoma cells with G9a knockdown. a-Tubulin levels are shown as the loading control. 
doi:1 0.1 371/journal.pone.01 06962.g006 



[40,41] (Fig. 3C). These results showed that G9a inhibition in 
neuroblastoma cells completely blocked cell proliferation, arresting 
the cells in the Gl phase. Together, these results suggest that G9a 
plays an important role in neuroblastoma cell growth and 
proliferation. 

Inhibition of G9a induces neuroblastoma cell autophagy 

To elucidate the contribution of G9a to the inhibition of 
neuroblastoma cell proliferation, we conducted immunofluores- 
cence staining. After 2 days of treatment, an autophagosome-like 
structure was observed in neuroblastoma cells, and this structure 
developed in a time-dependent manner after 4 days of treatment 



(Fig. 2 A). Next, we investigated LC3B expression, which is a 
marker of autophagy [42-44]. As shown in Fig. 4A and B, we 
found abundant LC3B expression in cells with G9a inhibition 
induced by BIX01294 treatment. To gain insight into the 
molecular mechanism underlying G9a inhibition-induced LC3B 
expression, we analyzed the protein expression of all five 
neuroblastoma cells. BIX01294 treatment had no significant effect 
on G9a expression in neuroblastoma cells, but it led to a marked 
down-regulation of H3K9me2, which was reported to be a major 
product in G9a-mediated H3K9 methylation [6] (Fig. 40). These 
findings indicate that the inhibitor BIX01294 had no effect on 
G9a protein expression, but it suppressed the methyltransferase 
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Figure 7. Downregulation of G9a induces autophagy and decreases tumorigenicity in neuroblastoma cells. A, immunofluorescence 
analysis of three neuroblastoma cells expressing GFPsi, G9asi#3 or G9asi#4. Scale bar, 5 |jm. B, statistical analysis of LC3B puncta in neuroblastoma 
cells expressing GFPsi, G9asi#3 or G9asi#4. Each column represents the average obtained from three independent experiments; error bars, SD. 
Statistical analysis was performed using two-tailed student's t-test, *p£0.01 . C and D, western blot analysis of G9a function (C) and autophagy-related 
genes (D) in neuroblastoma cells expressing GFPsi or G9asi#4. a-Tubulin levels are shown as the loading control. Cells expressing GFPsi are shown as 
the biological control. E, neuroblastoma cells were plated at 1x10^ cells per well in six-well culture plates. After 14 to 21 days of culture, soft agar 
colonies developed with cells expressing GFPsi. As shown, the cells with G9a knockdown were observed to give rise to small and scanty colonies in 
soft agar. Scale bars, 50 jim. F, Colonies that were larger than 0.5 mm or that contained more than 50 cells were recorded. Each column represents 
the average obtained from three independent experiments; error bars, SD. Statistical analysis was performed using two-tailed student's t-test, *p£ 
0.01. 

doi:1 0.1 371 /journal.pone.01 06962.g007 



function of G9a. Moreover, we investigated whether BIX01294 
treatment induces the expression of genes known to participate in 
autophagosome formation. The Western blot assay revealed tliat 
the expression of ATGs and LC3B was markedly up-regulated 
after treatment (Fig. 4D). These results indicate that the loss of 
G9a function induces autophagy and autophagosome formation in 
neuroblastoma cells. 

Inhibition of G9a decreases tumorigenicity of 
neuroblastoma cells 

To evaluate the role of G9a for the tumorigenicity of 
neuroblastoma cells, we first examined this phenomenon using 
soft agar clonogenic assays. The BE(2)-C cells were plated at 
1x10"* cells per well (six-well culture plates). After plating, 
immediate examination under a microscope revealed mosdy 
individual cells and colonies, defined as a collection of more than 
50 cells appearing after 14 to 21 days. The BE(2)-C cells treated 



with 5 |aM BIX01294 were observed to give rise to small and 
scanty colonies in soft agar compared to the cells treated with 
water. Similar results were also obtained with several other 
neuroblastoma cell lines, including SK-N-AS, SK-N-DZ, and SK- 
N-Fl cells (Fig. 5A and B). Together, these results indicate that 
G9a inhibition prevents the clonogenic activity of neuroblastoma 
cells. Next, we examined the effect of G9a on the ability of 
neuroblastoma cells to induce tumors in immunodeficient mice. 
The neuroblastoma cell lines BE(2)-C and SK-N-AS were injected 
subcutaneously into the flanks of NOD/SCID mice. After 1 week 
of tumor growth, BIX01294 was injected intraperitoneally at 
4 mg/kg. NOD/SCID mice were injected with water as a control. 
The mice without BIX01294 treatment developed large tumor 
masses after three weeks, whereas BIX01294 injection significandy 
diminished the tumorigenic activity of neuroblastoma cells during 
the same time period (Fig. 5G and D). AH of the above results 
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demonstrated that G9a inhibition decreases the tumorigenicity of 
neuroblastoma cells. 

Downregulation of G9a represses cell proliferation and 
induces autophagy in neuroblastoma 

In addition to the inhibition of G9a by BIX01294, we knocked 
down G9a in the three neuroblastoma cell lines SK-N-AS, BE(2)- 
C and SHEP 1 . We employed five lentiviral constructs expressing 
short hairpin RNA sequences against human G9a (G9asi), and 
three of them significantly reduced G9a expression in the 
neuroblastoma cells (Figure 6A). We performed the subsequent 
experiments with neuroblastoma cells expressing either G9asi#3 
or G9asi#4. Biochemically, we noted that G9a knockdown in the 
three neuroblastoma cell lines induced cell growth arrest, which 
was determined by morphology and proliferation assays (Fig. SIC, 
Fig. 6C and D). Western blot analysis showed that G9a 
knockdown also led to a marked down-regulation of CyclinDl, 
CDK4 and CDK6 in the three neuroblastoma cell lines (Fig. 6E, 
Fig. S2A). These results showed that G9a down-regulation 
represses neuroblastoma cell growth and proliferation. Further- 
more, we noted that G9a knockdown in three neuroblastoma cell 
lines also induced autophagosome formation and LC3B expression 
as determined by morphology and immunofluorescence assays 
(Fig. 6B, Fig. 7A and B). In addition, G9a knockdown in the 
neuroblastoma cells markedly inhibited the H3K9me2 level and 
up-regulated the ATGs and LC3B expression levels (Fig. 7C and 
D, Fig. S2B and C). These findings demonstrate that the down- 
regulation of G9a-induced autophagy in neuroblastoma cells. 
Moreover, we examined the tumorigenicity of the G9a-knock- 
down cells. The result showed that the neuroblastoma cells that 
underwent G9a knockdown failed to form large colonies in soft 
agar (Fig. 7E and F, Fig. S2D and E). Collectively, these results 
indicate that the knockdown of G9a in neuroblastoma cells 
represses cell proliferation, decreases tumorigenicity, and induces 
autophagy. 

Discussion 

The intent of this paper was to characterize the role of the 
H3K9 methyltransferase G9a in neuroblastoma cell growth and 
tumor development. In this study, we investigated the effect of G9a 
together with its novel small molecular inhibitor BIX01294 in 
neuroblastoma cells. Specifically, we sought to determine the 
function of G9a in cell growth, tumorigenicity and histone 
methylation in neuroblastoma cells. Our results demonstrated 
G9a expression in all five neuroblastoma cell lines, SK-N-AS, 
BE(2)-C, SK-N-DZ, SK-N-Fl, and SHEPl, suggesting a common 
expression of G9a in neuroblastoma. Furthermore, our results, 
based on two independent microarray databases, indicated that 
higher levels of G9a predict poor neuroblastoma prognosis. In 
addition, we found that tumors from older patients (>18 months) 
expressed G9a mRNA at significantly higher levels than tumors 
from younger patients (<18 months), suggesting that neuroblas- 
tomas initiated during early embryonic development are more 
likely to express lower levels of G9a. 

In addition, our data indicate that the treatment of neuroblas- 
toma cells with BIX01294 significantly inhibited cell growth and 
proliferation. The results from this study further demonstrated that 
BIX01294 treatment suppressed the tumorigenicity of neuroblas- 
toma cells both in vitro and in vivo. It was unclear how G9a 
inhibition caused by this small molecular agent promotes growth 
suppression through cell death or the inhibition of cell division 
[45,46]. To determine the mechanism by which the growth- 
inhibitory effects of G9a are overcome, flow cytometry experi- 



ments were performed. The results demonstrated that the cell 
cycle was altered after drug treatment. For example, the number 
of cells in Gl phase was increased and the number of cells in S 
phase was decreased. Our findings based on BIX01294-induced 
inhibition of G9a suggest that G9a may have an oncogenic 
function in the pathogenesis of neuroblastoma. 

Interestingly, we found several specific characteristics of 
autophagy after BIX01294 treatment, including the appearance 
of membranous vacuoles and the accumulation of LC3B in 
autophagosome [47-49]. Similar results were observed in G9a- 
knockdown cells. Recently, it was reported that BIX01294 induces 
apoptosis in three neuroblastoma cell lines (LAl-55n, IMR-5, and 
NMB) [50], which was different from our findings. Thus, it was 
particularly interesting to identify the molecular mechanism of the 
regulation mediated by G9a in neuroblastoma growth. BIX01294 
may induce autophagy-associated cell death in MCF-7 cells [51]; 
therefore, we speculate that G9a is involved in a novel regulatory 
pathway of autophagy in neuroblastoma. Herein, we provide 
evidence that the pharmacologic inhibition or genetic depletion of 
the histone methyltransferase G9a leads to the suppression of cell 
proliferation, the formation of autophagosome-hke structures, the 
aggregation of LC3B and the expression of ATGs. Moreover, our 
data suggest that G9a represses the level of histone H3K9me2, 
leading to increased expression of LC 3B and ATGs to support the 
process of autophagy. Taken together, these data indicate a role 
for G9a in the epigenetic regulation of autophagy and suggest that 
either pharmacologic manipulation of G9a or the genetic 
regulation of G9a may impact the autophagy- process. 

In conclusion, our results define the critical function of G9a in 
cell growth and tumorigenicity and suggest that G9a is involved in 
a novel regulatory pathway of autophagy in neuroblastoma. 

Supporting information 

Figure SI Inhibition of G9a suppresses neuroblastoma 
cell proliferation. A, neuroblastoma BE(2)-C cells were treated 
with BIX()1294 for 48 h at concentrations of 1, 2.5, 5, 7.5, and 
10 |tM. Proliferation was determined by trypan blue exclusion to 
differentiate between dead and live cells. Error bars, SD, n = 5. 
Statistical analysis was performed using two-tailed student's t-test, 
*pS0.01. B, morphologic examination of five neuroblastoma cells 
treated with 5 |iM BIX()1294 or water for 2 and 4 days. Scale 
bars, 20 |J,m. C, morphological examination of three neuroblas- 
toma cell lines express GFPsi, G9asi#3 or G9asi^4. Cells 
expressing GFPsi are shown as the biological control. Scale bar, 
20 |.lm. 
(TIF) 

Figure S2 Downregulation of G9a represses neuroblas- 
toma cell proliferation and tumorigenicity and induces 
autophagy. A, B and C, western blot analysis of cyclins and 
CDKs associated with the Gl phase (A), G9a function (B) and 
autophagy-related genes (C) in neuroblastoma cells expressing 
GFPsi or G9asi#3. ot-Tubulin levels are shown as the loading 
control. Cells expressing GFPsi are shown as the biological 
control. D, neuroblastoma cells were plated at 1 x 10^ cells per well 
in six-wen culture plates. After 14 to 21 days of culture, soft agar 
colonies grown with cells expressing GFPsi. As shown, the [:ells 
with G9a knockdown were observed to give rise to small and 
scanty colonies in soft agar. Scale bars, 50 |im. E, colonies that 
were larger than 0.5 mm or that contained more than 50 cells 
were recorded. Each column represents the average obtained from 
three independent experiments; error bars, SD. Statistical analysis 
was performed using two-tailed student's t-test, *p£0.01. 
(TIF) 
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